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Abstract 

We study the spin-dependent quark distributions at large transverse momentum. We derive their 
transverse momentum behaviors in the colhnear factorization approach in this region. We further 
calculate the angular distribution of the Drell-Yan lepton pair production with polarized beams 
and present the results in terms of the collinear twist-three quark- gluon correlation functions. In 
the intermediate transverse momentum region, we find that the two approaches: the collinear 
factorization and the transverse momentum dependent factorization approaches are consistent in 
the description of the lepton pair angular distributions. 
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I. INTRODUCTION 



Spin dependent semi-inclusive hadronic processes have attracted much interest from both 
experiment and theory sides in recent years. These processes provide us more opportunities 
to study the Quantum Chromodynamics (QCD) and internal structure of the hadrons, as 
compared to the inclusive hadronic processes or spin averaged processes. Measurements 
have been made in different reactions. In particular, the single transverse spin asymmetry 
(SSA) phenomena observed in various hadronic processes [l|, B , P H R| h ave stimulated 
remarkable theoretical developments 0, i, M, Cl, E, HM Elirlli, 13, M, M, H H 
23, 24, 2^, 2^. Among the theoretical developments, two approaches in the QCD framework 



have been most explored: the higher twist coUinear factorization approach 



and the transverse momentum dependent (TMD) approach [ll|, ll2| 
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29|,|30| 



s, the spin-dependent differential 
cross sections can be calculated in terms of the collinear twist three quark-gluon correlation 
functions in the collinear factorization formalism and the TMD distributions in the TMD 
factorization approach. Such functions generalize the original Feynman parton picture, 
where the partons only carry longitudinal momentum fraction of the parent hadron. They 
will certainly provide more information on hadron structure. Further study has shown that 
the transverse momentum dependence of the naive-time-reversal-odd TMD distributions 
which are responsible for the SSAs can be calculated and expressed in terms of the collinear 
twist-three correlation functions at the large transverse momentum. By using these results, 
it was shown that the above two approaches are consistent in the intermediate transverse 

In this paper, we will extend these studies to 



momentum region where both apply [31|, |32|, [33 



more general Drell-Yan processes, in particular the lepton angular distributions in polarized 
nucleon-nucleon scattering. 

The single transverse spin asymmetry in the Drell-Yan lepton pair production process has 



been used as an example to demonstrate the consistency between these two approaches [31 



where the single transverse spin asymmetry is represented as a correlation between the 
lepton pair transverse momentum q± and the transverse polarization vector S±. For this 
contribution, the transverse spin dependent differential cross section is proportional to 
da{S±) oc e^^ S\q'^. The transverse momentum of the lepton pair is also the transverse 
momentum of the virtual photon which decays into the lepton pair in the Drell-Yan process. 
Therefore, this SSA is related to the quark Sivers function, and is the only contribution from 
the quark-antiquark channel. However, if we further study the lepton angular distribution 
in the Drell-Yan process [13], it will open more contributions to the single spin asymmetries, 
as well as other spin dependent observables [35, 3^. More recently, by analyzing the general 
Lorentz structure of the hadronic tensor, the complete spin and transverse momentum de- 
pendent angular distribution of lepton pair has been presented in Ref . [22| , and 48 structure 
functions will contribute. One can calculate some of them which are leading power in the 



22 



context of the TMD factorization and relate them to the TMD distributions 

In this paper, we will study the angular distribution of the lepton pair in the polarized 
Drell-Yan process at large transverse momentum of the lepton pair. The relevant calculations 
are carried out in the collinear factorization framework. We mainly focus on the single spin 
asymmetry Aut and double spin asymmetry Alt-, taking into account the contributions from 
the twist-three quark-gluon correlation functions from one of the incident hadrons. We will 
limit ourselves up to this order in the calculations. For the single spin asymmetr y A tjt, the 



corresponding twist-three quark-gluon correlation functions are those studied in [28|, |29|, |3C 
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31 . 3^ and the calculations will be similar. On the other hand, for the Alt asymmetry, more 
rluon correlation functions will contribute and the calculations will be different 
m H, H HH, [13]. Following the same procedure as that in Refs. j3ll . 



general quark- 
from those in 

we will compare the predictions from the two formalisms and check their consistency. 

To pursue this aim, we will calculate the TMD quark distributions at large transverse 
momentum k±_ ^ Aqcd, and express them in terms of the collinear correlation func- 



tions. There are eight leading order TMD quark distributions three k±-even TMD 

distributions fi{x, k±),giL{x, k±), hi{x, k±); two naive-time-reversal-odd TMD quark distri- 
butions /fir(x, fcj,) (Sivers function), hj^{x,k±) (Boer-Mulders function); and three naive- 
time-reversal-even but fc_L-odd TMD quark distributions giT{x,k±), hiL{x, 
The transverse momentum dependence of these TMD quark distributions can be calculated 
from perturbative QCD in the collinear factorization framework. For the fc^-even TMD 
quark distributions, the results are well-known, and can be expressed in terms of the in- 
tegrated leading-twist parton distributions (see, e.g., |23|). The naive-time- reversal-odd 



isee, e.g., 

TMD quark distributions (the Sivers function and Boer-Mulders function) have also been 
calculated |3l|, 3^. In this paper, we will extend these calculations to the two naive-time- 
reversal-even but fc_L-odd TMD quark distributions gir andJiiL. These results will depend 
on the the novel twist three distributions g{x), h{x) (2"" 



and the general twist-three 



quark-gluon correlation functions Gd, Gd, Hd, and Ed [33]. The last TMD quark distribu- 



tion hj-rp{x, k 



will involve twist-four quark-gluon correlation functions. We will not discuss 
it in this paper. 

The rest of paper is organized as follows. In Sec. II, we give a brief review on the twist- 
three quark-gluon collinear correlation matrix elements and discuss the relation between 
them and the TMD distributions. General feature of the TMD quark distributions at large 
transverse momentum will be presented in Sec. III. In Sec. IV, we derive the naive-time- 
reversal-even TMD quark distributions girix, k±), hnix, k±) in the twist-three quark-gluon 
correlation approach. In Sec. V, we calculate the relevant polarized Drell-Yan differential 
cross section using the same collinear factorization and compare to the results from the 
TMD factorization. We conclude the paper in Sec. VI. 



II. TWIST-3 CORRELATION MATRIX ELEMENTS AND TMD DISTRIBU- 
TIONS 

In order to extract more information on hadron structure, various spin dependent and/or 
transverse momentum dependent parton correlation functions have been introduced based 
on the QCD factorization theorem. They are universal between SIDIS and Drell-Yan pro- 
cesses (up to a sign for the naive-time-reversal-odd TMD parton distributions) and can be 
pinned down by a complete set of experiments. It has been shown that there exits interest- 
ing connections between the twist-three collinear functions and the naive-time-reversal-odd 
TMD quark distributions j20l . 30| . In this section, we will review the general property of the 



collinear correlation functions and introduce two novel twist-three functions g,h [2(|. We 
will further explore their relations to the naive-time-reversal-even but k±-odd TMD quark 
distributions giT and hiL. 

Let us start by introducing the following collinear quark- ant iquark correlation matrix: 

Ma/^W = / ^e-^-''^y-{P,S\My')MOmS), (1) 
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where P, S are the hadron momentum and spin, respectively, and we have suppressed the 
hght-cone gauge hnks between different fields. The hadron momentum is proportional to 
the light cone vector = (l"*", 0~, 0^), whose conjugate light-cone vector is n = (0"*", 1~, 0^). 
X is the momentum fraction of the hadron carried by the quark. Up to twist-three level, the 
above matrix can be expanded as [s^], 



M{x) 



+ 



M 
2P+ 



e{x)l + gT{xh5^± + hL{x)-'y5[^, r/] 



(2) 



where A represents the helicity for the nucleon for the longitudinal polarized nucleon, S± is 
the transverse polarization vector, and M the hadron mass. The first three are the leading- 
twist quark distributions: spin average /i, longitudinal spin gi, and quark transversity hi 
distributions. The twist-three quark distributions: e(x), gri^), and do not have simple 

interpretations, and belong to more general quark-gluon correlation functions 37|. These 
correlation functions can be defined through the following matrix |37l . l38l . Hi, 



dy dVl -ixP+y- i(xi-x)P+yT 



27T 2tt 



{P,S\My-)zD^^{y^)MO)\P,S) , (3) 



where we have adopted the covariant derivative as iD'^ = id^ + gA''j_. The expansion of the 
above matrix contains the following four twist-three quark-gluon correlation functions. 



M^(a;,xi; 



+ 



M 

2P+ 
M 

2P+ 



Gd{x, xi)ie'^S±u^ + Gd{x, xi)S'^75|^ 
[Hd{x, xi)A757^|^ + Ed{x, xi)'j'^^] . 



(4) 



By imposing the hermiticity, parity and time-reversal invariance, we will have the following 
constrains. 



Gd{x, xi) = -Gd{xi,x), Gd{x, xi) = Gd{xi, x), 
Hd{x,xi) = Hd{xi,x), Ed{x,xi) = -Ed{xi,x) 



(5) 



and these functions are real. As mentioned, the twist-three quark distribution qrix), e(x), 



and hi{x) can be expressed in terms of the above quark-gluon correlation functions [37|, |41 



9t[x) 



X 



dxi 



Gd{x,Xi) + Gd{x,Xi 



hr.(x) = - 



dxiHi:){x, Xi) 
2 f 

e(x) = — dxiEoiXyXi) . 
X J 



(6) 
(7) 
(8) 



Therefore, Gd, Gd, Hd-, and Ej:, functions are more fundamental, which becomes evident 
when we study the scale evolution for the twist-three quark distributions, and the next-to- 
leading order perturbative corrections to the relevant cross sections 41. 42 . 
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In terms of the twist expansion, of course, D-type correlations are not the only ones at 
the twist-three level. One can also define a set of the F-type twist-3 correlation matrix 
elements, 



dy dy^ 
27r 27r 



-ixP+y i{xi-x)P+y^ 



{P,S\ij^{y~)gF^^{y^)ij^{Q)\P,S) . (9) 



Again, the expansion of the above matrix defines the following F-type quark-gluon correla- 
tion functions, 



+ 



M 

Y 

M 



(10) 



where, for convenience, we have used different normalization factors for Tp(x,xi) and 
T^\x,xi) as compared to Ref. jsi]. Hi, with a relative factor 2tiM. Similarly, the par- 
ity and time-reversal invariance implies. 



Tp(x, xi) = Tp(xi,x), Tf{x,xi) = — Tp(xi,x), 
fP{x,Xi) = -fl^\xi,x), Tjf^(a;,xi) =Tjf^(xi,a;) 



fill 



The F-type correlation functions are usually regarded as an alternative but not independent 
functions in the calculations to the inclusive DIS structure functions, such as qt structure 
function 421]. On the other hand, it has been found that the F-type correlation functions are 
more relevant for the single transverse spin asymmetry, and have been intensively studied 28|, 



29l . |30| ■ By using the equation of motion, these two types of the correlation functions can 



be related to each other [20|, |30 



Goix^xi) 
Gd{x,xi) 
Ed{x,Xi) 



P 
P 
P 



Hd{x,Xi) = P 



X 


— Xi 




1 


X 


— Xi 




1 , 


X 


— Xi 




1 



Tf{x,xi), 

Tf{x, Xi) + S{x — xi)g{x), 



T, 



Tp\x, Xi) + 6{x — Xi)h{x). 



where P stands for the principal value, and g, h are given by j2o|. 



(12) 
(13) 
(14) 
(15) 



3^, 



-ixP^y 



271 



{p,smy- 



dc-F^^{c))m\p,s) 



M 



g{x)S^-f5:^ + /i(x)A757^|^ 



(16) 



From the above results, we find that indeed the F-type and D-type correlation functions are 
not completely independent, and they form an over-complete set of functions. However, we 
still need g and h to completely describe the associated physics at this order, in particular, 
for the calculation performed in this paper. In the real calculations, we can either use D-type 
or F-type plus g and h as a. complete set of twist-three functions. 
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In the following, we will further reveal the physical meaning of g and h, and build the 
connection between them and the transverse momentum dependent quark distributions. The 
TMD parton distributions are important generalization of the conventional Feynman parton 
distributions. Because of additional dependence on the transverse momentum of partons, 
these distributions open more opportunities to study the partonic structure in nucleon. 
The nontrivial correlations between the parton transverse momentum and the polarization 
vectors of the parent nucleon or the quark itself provide novel consequence in the transverse 
component in the hadronic processes, for example, the single transverse spin asymmetry. 
Of course, upon integral over transverse momentum, these TMD parton distributions will 
naturally connect to the leading-twist and higher-twist parton distributions. In this paper, 
we will focus on the TMD quark distributions, which are relevant to the Drell-Yan lepton pair 
production. The TMD quark distributions can be defined through the following matrix 18|, 

mM, 



Ma/six, kj 



dy d^y± i^p+.y- 
(27r)3 



(17) 

where x is the longitudinal momentum fraction and fc^ the transverse momentum carried 
by the quark. The gauge link is along the direction represented by v which is conjugated 
to p. In the case that we need to regulate the light-cone singularities, we will use an off- 



light-cone vector: v 3> f ^ and v±_ = 0, and further define ( 



P) /v . Compared 



with the integrated parton distributions definition in the above, we find that the two quark 
fields are not only separated by light-cone distance but also by the transverse distance 
^x, which is conjugate to the transverse momentum of the quark k± 40]. Because of this 



difference, the additional transverse gauge link for the TMD parton distributions has to be 
contained in order to make the above quark matrix gauge invariant [l9|, and the gauge link 



direction depends on the process Il8|,ll9(. Since we will study the Drell-Yan lepton pair 



production in this paper, in the following we will adopt the TMD definition for this process 



and the gauge link will go to — oo [18|, ll9|]. The gauge link plays an essential role in the 
naive time-reversal-odd TMD parton distributions. 

The leading order expansion of the quark distribution matrix M. contains eight quark 



distributions 14, 15, 21 



M 



fi{x, ki_) i) + k±)a'"'kf,p„ + giL{x, A;±)A75 ^ 



1 - - 1 

+j^giT{x, /c_l)75 ¥^{kj_ ■ S_l) + —hiLXia^ulbP^k"^ + hi{x, ki_)ia^u'lbP^S''j_ 

+^hUx, k^)ta,,^,p^ (k^ ■ S^kl - hls{\ + l^fM^, k^)e^''"Pi,p,k^Sp 



(18) 



Out of the eight TMD distributions, three of them are associated with the A;_|_-even struc- 
ture: fi{x,k±), giL{x,k±), and hi{x,k±). They are simple extension of the above in- 
tegrated quark distributions. The other five distributions are associated with the k±- 
odd structures, and hence vanish when k± are integrated out for A^^/j. For an un- 
polarized nucleon target, one can introduce the unpolarized quark distribution fi{x,kj_) 
and naive-time-reversal-odd transversely-polarized quark distribution hj;{x,k±), the Boer- 
Mulders function. For a longitudinally-polarized nucleon, one introduces a longitudinally- 
polarized quark distribution giL{x, k±) and a transversely-polarized distribution /i^^(a;, k±). 
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Finally, for a transversely-polarized nucleon, one introduces a quark spin-independent dis- 
tribution f^{x,k±), the Sivers function, and a longitudinally-polarized quark polarization 
giT{x, k±), a symmetrical transversely-polarized quark distribution hi{x,k±) and an asym- 
metric transversely-polarized quark distribution h^rp[x, kx). 

If we weight the integral of M.ap with linear dependent transverse momentum, the k \ 



odd quark distributions will lead to the higher-twist quark-gluon correlation functions [20 
Four of them will correspond to the four quark-gluon correlation functions introduced above, 
including f^, h^, giT, hn. The last one hj-j,, as we mentioned, will correspond to a twist- 
four correlation function. First, the two naive-time-reversal-odd quark distributions and 
hj^ lead to the following quark-gluon correlations^ 



Tf{x,x) = I ^^/iilDY(x,fcx) , (19) 
Ti.^\x,x) = I ^^^hi\Mx,kx) , (20) 

where the TMD quark distributions follow their definitions in the Drell-Yan process. Sim- 
ilarly, the two naive-time-reversal-even but k±-odd quark distributions giT and hi^ can be 
related to the following twist-three matrix element, 

g[x) = / d'^k±^^giT{x,k±) , (21) 

h{x) = I £k^^^h,L{xM) . (22) 

Here, because they are naive-time-reversal-even distributions, they will not change sign 
between DIS and Drell-Yan processes. In summary, the four A;x-odd TMD quark distribu- 
tions correspond to the four twist-three quark-gluon correlation functions introduced above 
Eqs. (12)-(15). From these relations, we can further study the scale evolutions for the above 
twist three correlations 44, 



The above relations are only one aspect of the connections between the TMD quark 
distributions and higher-twist quark-gluon correlation functions. In the following section, we 
will explore another aspect, i.e., the large transverse momentum behavior for the TMD quark 
distributions in terms of the coUinear leading-twist or twist-three quark-gluon correlation 
functions. 



III. QUARK DISTRIBUTIONS AT LARGE TRANSVERSE MOMENTUM 

When the k^ is of the order of Aqqd, the TMD parton distribution functions are entirely 
non-perturbative objects. However, the transverse momentum dependence can be calculated 
in the perturbative QCD and related to the coUinear matrix elements as long as the k± is 
much larger than Aqcd, the typical nonperturbative scale. The collinear matrix elements 
are the relevant integrated leading-twist parton distributions or higher-twist quark-gluon 



These relations Eqs. (19-22) are valid at the leading order in perturbative expansion which we will use in 
this paper. They may differ from these forms at higher orders. 
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correlation functions. For example, for the three fc^-even quark distributions, they will de- 
pend on the integrated leading-twist parton distributions. However, for the fc_L-odd quark 
distributions, they depend on the twist-three (or twist-four) quark-gluon correlation func- 
tions. In general, we will have the following expression for the quark distributions at large 
transverse momentum 1231^. 



x' 



(23) 



where q{x, k±) represents the TMD quark distribution we are interested, fi represents the in- 
tegrated quark distribution for the fc^-even TMDs, and higher twist quark-gluon correlation 
function for the k±-odd TMDs. For the latter case, x' should be understood as two variables 
for the twist-three quark-gluon correlation functions as we discussed in the last section. The 
overall power behavior l/(/c^)" can be analyzed by the power counting rule The hard 
coefficient T-Cq/i{x; x') is calculated from perturbative QCD. In this paper, we will show the 
one-gluon radiation contribution to this hard coefficient. 

The fcj_-even TMD quark distribution functions, fi{x,k±), j,{x,k \ ), and hi{x,kj_) can 



be calculated from the associated integrated quark distributions 
contributions, they are expressed as 23l |. 



23|^. For the non-singlet 



fi{xB,k±) 
giiixB, kA_) 
hi{xB, k±) 



as 1 

as 1 
27c'kl 

as 1 



dx 

X 

dx 

X 

dx 

X 



fi{x) 



-9lL[X) 



fi{x) 



+ 5(l-0|ln^-l 



(24) 
(25) 
(26) 



where the color factor Cp = {N^-l)/2N^ with A^^ = 3, ^ = xb/x and = {2vPy/v^. Here, 
we have adopted an off-light-cone vector v to regulate the light-cone singularity associated 
with the above calculations [23| . 

In the same spirit, the naive-time- reversal-odd TMD distributions, the quark Sivers func- 
tion and Boer-Mulders function hj^ at large k± can be calculated perturbatively. The 
contributions come from the twist-three correlation matrix elements Tp, Tp, and Fur- 
thermore, it is known that the time-reversal odd TMD distributions is process dependent 
because the difference on the gauge link directions will lead to a sign difference between the 
SIDIS and Drell-Yan processes 17, 18, 19| . 



itIdy 



itIdis, 



hi 



DY 



-hi 



DIS 



(27) 



^ This is not a rigorous factorization formula. However, we shall be able to construct a QCD factorization 
formalism in the impact parameter 5-space for the TMD distributions 



23 



4C 



47| 



^ Mixing with the gluonic contributions will have to be taken into account for /i and gi distributions. In 
this paper, we will not discuss these contributions. 



8 



The quark Sivers function and Boer- Mulders function have been calculated 31, 



flT\'CiY{.XBj k±) 

hi\DYixB,k±) 
where the A factors are defined as 

1+e c 



as 


M2 


f dx 


vr 


{kiy J 


X 


as 




r dx 


IT 




X 



A,.+CpTF{x,x)5il-0\^n- 




A 



A 



hi 



CfTf{x, x) 



CfTP{x,x)- 



-0- 

2e 



0- 



_A 

2 

Ca 
2 



2 



Tf(x,xb) 



^-^T^^\x,xb) - ■^^^T^^'{x,x) 



2e 



Tf(x, x) 



TF{xB,x\?fQ) 

(31) 



where the color-factor Ca = N^.- 

From the above results, we can see that the large transverse momentum TMD quark 
distributions have a generic structure. They contain two parts: one part is similar to the 
splitting kernel for the relevant collinear functions, and one term is a delta function associated 
with a large logarithm hiC'^/k']_ which comes from the light-cone singularity regulated by an 
off-light cone gauge link discussed above. The splitting kernel may be different for different 
TMD quark distributions. However, the logarithmic term is the same for all of them. This 
is because this term is related to the soft gluon radiation and is spin- independent. We can 
also use this as an important consistent check for all the calculations. In the next section, 
we will calculate the two k^-odd but time-reversal even quark distributions Qit and 
and we will find that they have the same structure. 

The energy dependence of the TMD quark distributions, the derivative respected to C^, 
is controlled by the so-called CoUins-Soper evolution equation |40l . Si- These evolution 
equations can be used to perform soft-gluon resummation for the final k±_ distribution in 
the cross section and the TMD quark distributions. It is more convenient to study this 
resummation in the impact parameter 6-space liT 
future. 



49 . We will address this issue in the 



IV. TRANSVERSE MOMENTUM DEPENDENT QUARK DISTRIBUTIONS qit 

AND hiL 

In this section, we will calculate the large transverse momentum behavior for the two 
naive-time-reversal-even but fc^-odd quark distributions: Qit and hn. These calculations 
will follow the previous calculations on the naive-time-reversal-odd quark distributions /fi^ 
and h^. However, they differ in a significant way. As shown in (itI . 18 , 19| . for the naive- 



time-reversal-odd distributions, the gauge link contributions play very important roles. For 
example, these time-reversal-odd distributions will have opposite signs between SIDIS and 



^ The derivative terms in these results |31| have been transformed into the non-derivative terms by partial 
integrals. The associated boundary terms were canceled out by the same boundary terms from the 
derivative terms 14511. 
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(a) 




FIG. 1: Generic diagram interpretations for the twist expansions in the high energy scattering am- 
plitudes up to twist-three level: (a) corresponds to a leading twist matrix element {ipip); (h)-(d) for 
twist-three contributions, (h) for {tpd±ip), (c) for {tpA^ip), and (d) for {ipd^A^ip) . Additional 
gluon connection between hard partonic part and the non-perturbative nucleon structure part can be 
added to these diagrams. This is because they do not change the power counting in these diagrams. 
The contributions from these diagrams (b-d) are not gauge invariant individually. However, they 
will combine into the gauge invariant results in terms of the correlation functions introduced in 
Sec. II. 



Drell-Yan processes, because the gauge links go different directions. In the practical cal- 
culations, we have to take pole contributions from the gauge links in these TMD quark 
distributions 3l|, |32j. The calculations for qit and hiL are different. Because they are 
naive-time-reversal-even, we do not take the pole contributions from the gauge links. That 
makes the calculations a little more involved, as the two- variable dependent correlation func- 
tions will enter explicitly. As we mentioned above, these two TMD quark distributions will 
depend on the correlation functions. Go (Tf), Gd (Tf), and Hd {Tp'^). Moreover, they will 
have contributions from the twist-three function g and h. In the following, we will calculate 
these contributions. 

The twist expansion will be the key step in the calculations. The tech niq ue used to 
calculate these contributions has been well developed in the last few decades [28|, 



3i,[3i,|33,[3i,|4ll,|4i. In the following, we will sketch the main points of our calculations for 
the TMD quark distributions Qit and hi^ and the Drell-Yan lepton pair production cross 
section in the next section. 

In the twist expansion, a set of non-perturbative matrix elements of the hadron state 
will be analyzed according to the power counting of the associated contributions. At the 
twist-three order, from a generic power counting we have contributions from the following 
matrix elements [3a, [ 



42| 



(32) 



where the quark field spin indices have been suppressed for simplicity. The relevant Feynman 
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diagrams can be drawn accordingly. We illustrate the typical diagrams for the associated 
contributions from the above matrix elements in Fig. 1. For comparison, we have also shown 
the diagram corresponding to the leading-twist contribution from the matrix element {i^'ip) 
in Fig. la. Figs. Ib-d represent the contributions up to twist-three quark-gluon correlation 
matrix elements. Fig. lb corresponds to the contributions from the matrix element {ipd±ilj) , 
Fig. Ic from {ijjAj_ip), and Fig. Id from {ipd±A'^il)) . Because of additional gluon compo- 
nent in the matrix elements for Fig. Ic and d, there will be gluon attachment from the 
nonperturbative part to the perturbative part as shown in these diagrams. To calculate the 
contributions from Fig. lb and d, we have to do collinear expansion of the partonic scattering 
amplitudes in terms of and k^^ = pgj. ~ respectively. These expansions, combining 
with the quark field and gluon field, will lead to the contributions in terms of the matrix 
elements: {ipd±ijj), and {ipd^A^ip) . The calculation of Fig. lb is straightforward, without 
expansion in terms of the transverse momenta of the quarks and gluon. Furthermore, all 
these calculations have to be combined into the gauge invariant matrix elements, such as 
Cd, Gd; H£i, Ed, Tp, Tp, g, and h. 

However, these functions form an over-complete set of correlation functions at this order 
(twist-three) as we discussed in Sec. II. Therefore, we can express the results in terms of 
either D-tjpe or F-type correlation functions. For example, in the calculations of the twist- 



three gx structure function [4l|, |42|, one has to combine the contributions from {ipdjjp) 
and {ipA±ilj) into the gauge invariant form {ipD^ip) which is associated with Gp, and Gp. 
Meanwhile, for the single spin asymmetry observables (or the naive-time-reversal-odd TMD 



quark distributions) [28|, l29|, l30|, iMl, l32| , it is more convenient to calculate the contributions 
in terms of {il!d±A'^ilj) matrix element. Such matrix element is part of the gauge invariant 
matrix element {tpF^^ip) which is associated with the twist-three correlation functions Tp, 
Tp, and T^p\ The contributions from the diagrams associated with {il)d^ Ax_il)) have also 
been shown to exactly coincide with those from {iljd±A'^ip) to form a complete result into 
the form in terms of (ipF-^^ip) [30] . 

In the above two examples, it seems that the g and h functions are not necessary in 
these calculations, because they do not appear in the final results. This is only because 
in these calculations one has chosen a particular set of correlation functions for the final 
results. Otherwise, g and h functions will show up if we choose different set of correlation 
functions. For example, the structure function gp can be solely expressed in terms of 
and Gp. However, if we rewrite gp structure function in terms of Tp and Tp, we will have to 
introduce the g function, because of the relation of Eq. (13). Similar arguments apply to the 
single spin asymmetry calculations. Furthermore, the necessary of g and h will be manifest 
in the following calculations of the TMD quark distributions gip and hn. As shown below, 
their roles become so essential that we have to include them in the first place. This can also 
be seen from the relations between gip (Z^-il) and g [h) discussed in Sec. II. 

We will take gip calculation as an example to show how we perform the computation at 
twist-three level with the quark gluon correlation functions Gp {Tp), Gp {Tp), and g. The 
calculations for the TMD quark distribution hn and the Drell-Yan cross sections in the next 
section will follow accordingly. As outlined above, we first draw the associated Feynman 
diagrams for the large transverse momentum gip quark distribution. In order to calculate 
the large transverse momentum behavior for the gip function, we have to radiate a hard 
gluon. The relevant diagrams are plotted in Fig. 2, where the probing quark carries the 
momentum k = XpP + k± and the nucleon momentum is denoted by P. The double lines in 
these diagrams represent the gauge link expansion from the quark distribution definition in 
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(61) (62) (63) (64) 




(cl) (c2) (c3) (c4) 




(dl) (d2) (d3) (d4) 




(el) (e2) (e3) (e4) 



FIG. 2: Feynman diagrams for the TMD quark distributions at large transverse momentum cal- 
culated from the twist-three quark-gluon correlation functions. The mirror diagrams of (bl)-(e4) 
are not shown, but included in the final results. (al)-(a4) correspond to the contributions from 
the matrix elements of {'ipd±tp); (bl)-(b4), (el)-(e4), (cl) and (c3) correspond to the diagrams 
contributions from {'ipA±'ip); and (bl)-(e4) for {tpdj_A'^^) . 

Eq. (17). Again, these diagrams include the contributions from the matrix element {ipd±ip) 
(al-a4); from (ipA'^ip) (bl-e4). To obtain a complete result, we have to attach the gluon to 
all possible places as shown in the diagrams (bl-e4). This also guarantees that we will get 
the gauge invariant result. The mirror diagrams of (bl)-(e4) where the gluon attaches to 
the right of the cutting line of these diagrams were not shown in Fig. 2, but included in the 
final results. Part of the diagrams (bl)-(e4) correspond to the contributions from {ijjA±ilj), 
whereas all of them contribute to that from {iljd±A~^ilj). 

To calculate the TMD quark distribution qit at large transverse momentum, we first 
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compute the individual contributions from the matrix elements shown in Eq. (32). Then we 
will combine the individual results into the gauge invariant twist-three quark-gluon corre- 
lation functions defined in Sec. II. We can parameterize the associated matrix elements in 
Eq. (32) which correspond to our calculations. For example, the relevant {ipd^i/j) matrix 
elements are parameterized as 



Mi;(x) ^ 



dy- 



y'{p,smy-)id^^m\P,s) 



2tx 



(33) 



where Tq^ and Tq^ correspond to the parts in the gauge invariant functions Gd and G^, 
respectively. Similarly, we can define the associated {'il)A±_%p) matrix elements. 



dy dy^ 



271 2ti 



M 
2P+ 



Ta^ (x, xi)ie'^S±^^ + Ta^ (x, xi)Sl-f5^ 



(34) 



Notice that because of additional gluon attachment to the nucleon state, the above matrix 
elements will depend on two variables (x,Xi) which represent the momentum fractions car- 
ried by the quarks from left and right sides of cut line in the diagrams. In the case that 
there is no gluon attachment like that in Eq. (33), they are equal and become one variable. 
From the above matrix elements we can easily define those of {ipd^ A±_il)), 



M^+^ {x,xi) 



dy dy^ 



M 

Y 



27C 27r 

Tg+A^ (x, Xi)eYS±^^ + Tq+a^ (x, Xi)zS'^75j^ 

At the same order, we shall also have the following matrix elements, 

dy- dy^ 



(35) 



Mt^ . + {x,xi] 



2tt 2n 



-i^P+y- ei(^^-^)P^yi {P, S\^{y-)gd^^A+{y^)^{0)\P, S) 



M 



TQ^A+{x,xi)t^S^yi) + TQ^A+{x,xi)iS^^'^^'i) 



(36) 



From the definitions in Sec. II, we will find that the above matrix elements can form the 
following gauge invariant correlation functions. 



Tp(x,Xi) = Tq^a+{.X,Xi) +Tq+a^{x,xi) , 
Tf{x,Xi) = Td^A+ix,Xi) +TQ+A^ix,Xi) , 

g{x) = fo^{x)+ I dxiP — ^ — 



X — Xi 



d^A+[X,Xi, 



(37) 
(38) 

(39) 



where we have used the time-reversal invariance to derive the last equation. There is no 
similar relation for Tg^, which on the other hand can be related to Tf{x,x) depending on 
the choice of the boundary conditions for the gauge potential |50] , 



Adv : -Tp(x,x), A_l(+oo) = 
Tfl, (x) = ^ Ret : Tf{x,x), A^{-oo) = 

PV: 0, A^_{+oo) + A^{-oo) 



(40) 
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It has been shown that the final results on the single spin asymmetries will not depend 
on the boundary conditions for the gauge potential, although they correspond to different 
relations between the matrix element Tdj^{x) and Tp(x, x), and different contributions from 
individual diagrams H^l- 

Having sorted out the above relations, it is relative straightforward to perform the calcu- 
lations. As outlined above, we will calculate the Feynman diagram contributions in terms of 
the matrix elements at the right hand sides of Eqs. (37-39). These results will be combined 
into the gauge invariant correlation functions at the left sides of Eqs. (37-39). In the fol- 
lowing, we will calculate these contributions separately, and show that how we will combine 
them into the gauge invariant results. 



A. Contributions from Tq^^ and Tq^^ 

The contributions from Tg^ and Tg^ come from the diagrams Figs. 2(al-a4). As mentioned 
above, we will perform the collinear expansion to calculate their contributions. That is, the 
hard partonic part illustrated in the upper parts of these diagrams can be expanded in terms 
of the transverse momentum of the quark connecting to the nucleon state in these diagrams 
(see also Fig. lb). This momentum can be parameterized as 

P^ = xP^'+pI^ (41) 

where x is the longitudinal momentum fraction of the nucleon and p± is the transverse 
momentum. In the collinear expansion, we take the approximation that p_l <^ and only 
keep the leading non-trivial terms which are relevant for our calculations. For example, the 
hard partonic part H can be expanded as 

H{k,p)=H{k,p)\,=,P + pl^^^^\,=,P + --- , (42) 

where ellipsis stands for higher order expansion terms, a is a transverse index a = 1,2. 
The first term in the above expansion does not contribute to the TMD quark distribution 
OirixB, k_i_) at large transverse momentum. The second term will lead to the contribution 
from the matrix element Tg^ and Tg^. 

In the calculations, we substitute Eq. (42) into the hard partonic part H{k,p) in the 
Feynman diagrams of Figs. 2(al-a4), and take the linear term in the expansion. One par- 
ticular contribution is the so-called derivative term, which comes from the expansion of the 
on-shell condition for the radiated gluon ki = p — k. To calculate this contribution, we only 
keep the pj_ dependence in the delta function of the on-shell condition, and set p± = for all 
other factors in the hard partonic scattering amplitude. The final result will be proportional 
to the corresponding Born diagram in the collinear limit j2^ , 



gMxB,k^)\l = -,^Cp [ - f-x^fg^ix) ) (l + e) , (43) 



1^ 

-^Sx tt"^ kj_ J X \ dx 



where ^ = xb/x and 1/A;^ behavior comes from the power counting for the fc^-odd TMD 
quark distributions. 

For the non- derivative terms, we keep all p± dependence in the hard partonic scattering 
amplitude, and expand to the linear term in p±. Although it is tedious, the calculation is 
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straightforward, and we obtain 
1 



71 



(44) 



where ('^ has been introduced in Sec. III. After partial integrating for the derivative term, 
we can add these two terms Eqs. (43) and (44) together^. 



giT{xB,k^)\f^ 



7r2 



X 



(1-0- 



+ 5(1-0 In 



- 1 



(45) 



Similar calculations can be performed for the contributions from Tg^, and we find that it 
does not contribute to the TMD quark distribution giT^XB, k±). 



B. Contributions from Taj_ {Tq+a^ and T^^ {'^d+A±) 

Because the attaching gluon is transversely polarized {A±), the contributions from the 
matrix elements T^^ and T^^ come from the diagrams Figs. 2(bl-b4), (el-e4), (cl) and 
(c3). To calculate the contributions from these diagrams, we take the kinematics illustrated 
in Fig. 1(c), where the quark and gluon lines connecting the hard partonic part and the 
nucleon state only contain collinear momenta, 

K = XiP^ pt^ = X2P^, fc^^ = (x2-Xl)P^ (46) 

where kg is the attaching gluon momentum. These calculations are straightforward, and we 
obtain the contributions from Ta^ , 

as I f dxdxi- \ f ^ f , 2x| xb , 

giT{XB,k±)\f = —jj / Ta^X^Xi) <Cf\^ + 1 

-^-L -K"^ k\ J X [ \ Xi XiX X 

^ Ca (x| + xxi){2xb - X- xi) \ 
2 {xb — xi){x — Xi)xi J 

where we have used the symmetric property for Ta^x^xi) to combine the results from 
Fig. 2 and their mirrors. Because of additional gluon attachment in these diagrams, we will 
have two contributions from two different color factors, Cp and Ca- Similarly, we have the 
contribution from Ta^, 



giT{xB,k±)\TA^ = j dxdxiTAAx,xi) !^Cf 



2 

Xb _^ _ Xb _ ^ 

X2 If -I If 



2 (xi — xb)xi 



(48) 



Moreover, the above results can be translated into the contributions from Tq+a^ and Tq+a^- 
This is because we have the following relations between the above matrix elements, 

Ta^(x,xi)=P — ^ — fQ+AAx,xi), TaAx,xi) = P — - — T9+aAx,xi) , (49) 

X — Xi X — Xi 



^ Note that the boundary term when we partial integrate the derivative contribution was canceled out by 
the boundary term when we compute the derivative term in Eq. (43) [45i] . 
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where the imaginary parts in the right hand sides of the above equations have been dropped, 
because they do not contribute to the qit and hn calculations here. However, when we 
calculate the single spin asymmetry observables (such as the time- reversal-odd Sivers and 



Boer- Mulders functions), we have to keep these imaginary parts in the above equations [50 
Substituting the above results into Eqs. (47,48), we shall obtain the contributions from the 
matrix elements Tq+j^^ and Tq+j^^^. 

C. Contributions from Tg^A+ '^dj^A+ 

For these contributions, it is the A'^ component connecting from the nucleon state to the 
hard partonic part, and the gluon can attach to the gauge links in the Feynman diagrams. 
Therefore, we will have all diagrams in Figs. 2(bl-e4) contributing to the final results. More- 
over, since these matrix elements involve d±A'^, we have to perform the collinear expansion 
of the hard partonic parts in terms of the gluon transverse momentum. In doing so, we keep 
both transverse momenta for the two quark lines connecting the hard part and the nucleon 
state, 

pi; = xiP^ + p^^, = + p'i^ . (50) 

Clearly, the kinematics tell us that fc^ = (x2 — xi)P^ + and = ~ Pi±- The 
corresponding collinear expansion of the hard partonic part takes the following form. 



Hik;pi,p2) = H{k,pi,p2) 



IPlX=P2±=0 



dH{k;pi,p2) ^ dH{k;pi,p2) ^ 

Again, the expansion coefficients in the above equation can be calculated following the same 
method as we discussed in the above for that for the Tq^ contribution. For example, there 
is also derivative terms associated with Tq^a+ matrix element. This contribution also comes 
from the expansion of the delta function for the on-shell condition for the radiated gluon ki. 
The derivation for this part is similar, and we obtain the following result, 

gMxBM)\l^^^ = ^2^^C, I I dx,P^Te,M^,x,] 

(52) 

where the same hard coefficient as that for Tg appears as it should be due to the gauge 
invariance. As we mentioned above, this derivative term comes from the delta function 
expansion for the on-shell condition of ki. To calculate this contribution, we only keep the 
Pi± dependence in this delta function, and set pi± = for all other factors in the hard 
partonic amplitude. Because of this and the fact that it is the A'^ insertion in the hard part, 
we can use the Ward identity argument to summarize all diagrams into a simple factorization 
form: a product of hard partonic part in the Born diagram without the gluon insertion and 
the factor l/{x — Xi) representing the gluon insertion. 

Because they have the same hard coefficient, we can combine the derivative contributions 
from Tg and Tg^A+ together. In particular, by using Eq. (39), we can add the results from 
Eqs. (43) and (52), 
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in term of the gauge invariant correlation function g. This shows how we obtain the gauge 
invariant results from the individual contributions. Moreover, it also demonstrates that g is 
an independent contribution to the large transverse momentum gix quark distribution. 

The non-derivative contributions from Tq^a+ can be calculated accordingly, by keeping 
all the dependence in the partonic amplitude. The final result is. 



girixB, k_i 



7r2 fcl 



dxdxi 



X X — Xi 

Ca {xl + .r.ri){2.ru - .r - 



+ Cf 



X 



X^ X\ 



2x'|j X]^ ^ 



X\X 



X 



+ 



{xb - X]){x - X])Xx 



(54) 



In the above results, there are two terms with the color-factor Cp- Clearly, one term will 
be combined with that of Tq in Eq. (44) to form the contribution from the gauge invariant 
function g[x\ similar to the case for the above derivative contributions. The other term will 
be combined with that of Tg+Ax from Eqs. (47,49) to form the contribution from the gauge 
invariant function Tp{x^x\). 

Similarly, we can calculate the contributions from the matrix element Tq^^-^ , 



gvr{xB-,'k\)\T, 



a.. 



7r2 k\ 



dxdxi 1 \ \ ( , Xb 
-Td^A+{x, xi) <^Cf { — + — 



Xl 



x^ 



Xl 



X 



+ 



3^3^ 1 

2 (xi - xb)xi 



(55) 



Again, this result will combine with that from Tq+a^ from Eq. (48,49) to form the gauge 
invariant result in terms of the gauge invariant function Tf{x,Xi). 

Combining all these results together, we will obtain the final results for the TMD quark 
distribution giT{xB,k±) at large transverse momentum. 



giT{xB,k±) 



as M2 
^(fci)2 



dx 

X 



A,^^ + CFg{x)5{C-l) 



(56) 



where Ag^r^ is given by. 



-I' 



dxi \ g{x)CF 



X 



X 



B . 



Ix^ 



5{x-xi) +P 



X 



+p- 



X — Xl 



Xl XiX 

Tf{x,xi) 



Xb 
X 



1 + 



Tf{x,Xi) 

X — Xl 

Ca {x% + xxi){2xb — x — xi) 
2 {xb — xi){x — xi)xi 



( % + ^ _ ^ _ 1 



x^ 



Xl 



X 



Ca X xxi 
2 {xi - xb)xi 



(57) 



and we have partial integrated the derivative terms. Using the identities Eqs. (12) and (13), 
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Ag-^^ can be transformed into the D-type correlation functions, 



A 



giT 



dxi < g{x 



1+e CAi+e 



+ Gd{x,xi] 
+ Gd{x,xi) 



X 



— - 1 



5{x — xi) 

Ca {x\ + xxi){2xb — X — Xi) 



Ct. 



[xb — a;i)(x — x\)x\ 



C A X B XX\ 

2 (xi — 



(58) 



This result can also be obtained by directly combining the contributions from Tq^ in Eq. (45), 
Ta^ in Eq. (47), and Ta^ in Eq. (48). 

For TMD distribution hulx, k±), the perturbative calculation follow the similar proce- 
dure. It receives contributions from Tjj^\x, xi), h{x). We skip the detailed derivation, and 
list the final result, 



hiLixB, ki 



a., [ dx 

X 



2 (l.'i \2 



Ah,^+CFh{x)5{i-l)[\n 



where Ah^j^ is defined by. 



Ah, 



dxi { Cp 



h{x) 



2e 



4c! 



5(xi - x) + P-^f^/\x, xi] 



(59) 



X Tp (x, Xi , 



2(x — Xi — Xb) 2xb(xbx + XfiXi — x^ — x^) 

CiT' h 



Xi 2 (x^ — xi)(x — Xi)xi 

Similarly, it can be expressed in terms of the D-type function as, 

2e Ca 2i 



Ah, 



dxi < h{x) 



+ Hd{x,Xi) 



Cf 



Cp-, ^ S(xi — x) 

2(x — Xi — x^) 2xb(xbx + xbXi — x^ — x^) 



Xi 



+ 



{xb — xi)(x — Xi)xi 



(61) 



Indeed, the above results for the naive-time-reversal-even but k±-odd distributions also have 
the same structure as those discussed in Sec. III. 



V. SINGLE SPIN AuT AND DOUBLE SPIN Alt ASYMMETRIES IN THE 
DRELL-YAN LEPTON PAIR PRODUCTION PROCESS 

In this section, we will calculate the angular distribution of the lepton pair in the polarized 
Drell-Yan process, especially the Ajjt and Alt asymmetries: One of the incident hadrons is 
transverse polarized and another is unpolarized or longitudinal polarized . We focus on the 
lepton pair production in hadronic scattering, 

Ha + Hb^Y+X ^ tr + X , (62) 



which comes from the virtual photon decays. In the leading order, virtual photon is 



produced 



through quark-antiquark annihilation process, qq 7* in the parton picture [34 1. In the 
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rest frame of the lepton pair, we can define two angles [35| : one is the polar angle 9 between 
one lepton momentum and the incident hadron; the azimuthal angle is defined as the 
angle between the hadronic plane and the lepton plane. The general formalism has been 
worked out for the angular distributions in the polarized Drell-Yan process l^]. For our 
calculations of Aut and Alt-, one can write down the following general structure of the 
angular distribution of lepton pair in the Collins-Soper (CS) frame 2^ . 



da. 



{UT,LT) 



at. 



X 



d^qdn 2(27r)452Q2 
{\SaT\ [sin 0a ((1 + cos^ e)W^^ + (1 - cos^ e)WE'^ + sin 29 cos (j)W^'^ + sin^ 6 cos 2(f)W^l) 
+ cos 0a (sin 29 sin (j)W'f^ + sin^ 9 sin 20W^2a ) 
+ \SaT\SbL [cos 0a ((1 + cosH)W^^ + (1 - cos^ ^ + sm 29 COS (j)W + sin^9cos2(f)W^l) 

\ (63) 



+ sin 0a ( sin 29 sin (pW'^'^ + sin^ 9 sin 201^^"^^ 



'LT 



where the orientation of the transverse polarization of the hadron a is expressed through the 
CS-angle 0a, and and 9 have been introduced above. In the above expressions, SaT and 
SbL are the hadron a transverse polarization vector and hadron b longitudinal polarization 
vector, respectively. The angular-integrated cross section is expressed in terms of the W^^, 



WY^, W^^ and Wt^ as l22 



da. 



{UT,LT) 



a: 



127r352Q2 

{\SaT\ sin 0a {2WJY^ + W^^) + |5ar|5bLCOS0a {2W, 



(64) 



In particular, the structure function 2Wt'^ + Wi^^ has been calculated 3l|] , which represents 
the Sivers contribution to the single transverse spin asymmetry. From the above expression, 
we also see that the Sivers contribution is the only contribution to the single spin asymmetry 
when we integrate out the lepton angles. 

In the following, we will calculate the structure functions W^'^ and W^'^ in Eq. (63) 
in the intermediate transverse momentum region Aqcd (1± Q, and will compare the 
predictions from the collinear factorization and the transverse momentum dependent ap- 
proaches. For the single spin asymmetry W^'^, we follow the previous calculations [3l|. The 
only difference is that the calculation in 3l|] is equivalent to the virtual photon production 
in single transversely polarized nucleon-nucleon scattering, whereas we will contract the 
hadronic tensor to the lepton tensor to obtain the angular distributions of the lepton pair 
in this process. However, the technique method is the same. Again, as discussed in |31| . 
we will have soft pole and hard pole contributions ^, and they will have cancelation in the 
intermediated transverse momentum region. For the double spin asymmetry part W^^, we 
will follow the procedure in Sec. IV to calculate the twist-three contributions to the angular 
distributions. The corresponding correlation functions will be Gd (Tp), Gc (Tp), g, and so 
on. The similar Feynman diagrams can be drawn accordingly, which can be organized into 



^ The so-called soft-fermion pole will also contribute to the cross sections 3l|, |5l| which have been neglected 
in our calculations. However, these contributions do not change the conclusions of the consistency between 
the two approaches 31 1. 
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the contributions from the twist-three matrix elements {ipd±_'4)), {ipA^il)), and {ipd±A^ip) . 
These contributions are then combined into the contributions for the gauge invariant cor- 
relation functions Gd-, Gd, g, and so on. In order to calculate the different terms in the 
angular distribution Eq. (62), we choose to work in the CoUins-Soper frame [35|, where four 
orthogonal unit vectors are defined as [52|, |53| , 



■- 



Q 



Q± VQ^ + Ql 



qpP^ - qpP^ 



(65) 



where q^ is the virtual photon momentum, P, P are the momentum of two hadrons, and 
we further define = [P^ - (P ■ g)/g^g^]/\/5, = [P^ - (P ■ q)/q^q^]/VS with qp = 
P ■ q/ ^/S, qp = P ■ / a/S* and S the total hadron center of mass energy square, respectively. 
The structure function can be obtained by contracting the hadronic tensor W^"^ with six 



symmetric tensors constructed by the four orthogonal vectors [52, l53 



Z'^Z'^W, 



flU 1 



iZ^'X" + Z''X'')W^ 



^IV 1 



W^A = -{-X^^X" + Y^'Y'')W^, 



UT 



.{Y^'Z" + Y''Z^')W^^ , 



W'f^ = -(F^X*^ + Y''X^')W^ 



flU 5 



(66) 



and similar expressions hold for the UT and LT structure functions. 

Furthermore, we are interested in the cross section contributions in the intermediate 
transverse momentum region, Aqcd Q± ^ Q- To obtain the leading order contributions, 
we only keep the leading terms in Q±/Q, and neglect all higher order terms. With this power 
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counting expansion, six leading order structure functions survive and can be simplified as 



vrQi ^ Nc J X z 



dx dz 



A^.{z)5{l-i) + B^.{z)5{l-i)\h,{x) 



(67) 
(68) 



W, 



'UT 
AA 



/ii(a:) , (69) 



W, 



AA 



AA 



asM 



i(x) . (71) 



^2Q3 2^nJ 



dx dz 
X z 



A,,^{z)6{l-0 + B,,,{z)6{l-i)\h,{x) , (72) 



where ^ = xb/x, ^ = zb/z, xb and zb are defined as = Q/VSe^ and = QI^~Se~^ 
with y the rapidity of the lepton pair in the center of mass frame, respectively. The functions 
A/^j,, v4;j±, Agjy, Ah^^ have been defined in Sec. Ill with appropriate variable replacements. 
The functions 5y±,, 5^x, Bg^^, B^^^ are given by. 



i?,x(^) = CpT^^z.z) 
Bg^^(x) = Cpgix) 
Bh^^z) = CfHz) 



1+e 

(1-0+ 

.(1-0+ 

i+e^ 



2'5(e - l)ln^ 



- l)ln^ 



1-0- 

2^ 



L(i-0+ 



+ 25(e-l)ln^ 
+ 25(e - l)ln^ 



(73) 
(74) 
(75) 
(76) 



respectively. 

Meanwhile, the transverse momentum dependent factorization can be applied in the small 
transverse momentum, Q± <^ Q. The relevant structure functions can be written as 2^ . 



q± ■ k 



a± 



LT 



Ma 

g± • hi. 

Ma 
q± ■ ka± 

Ma 

q± ■ h± 



flT{xB,ka±)fl{zB,h±) , 

hi{xB, ka±)h^{zB, kb±) 



QirixB, ka±)giL{zB, kb±) 



hl{xB, ka±)hiL{zB, k 



b± 



(77) 
(78) 
(79) 
(80) 
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where the simple integral symbol represents a complicated integral 

J c ^ J 

and S{\i_), H{Q'^) are the soft factor and hard factor, respectively. In addition, the combina- 

tions of structure functions ^ , ^ also receive the leading power contribution 
from the product of TMD distributions /i^^ x and h^rp x h^j^ respectively, which is however 
beyond the scope of the present paper. 

When the transverse momenta kai_ and /cf,_L are of order Aqcd, the TMD distribution 
functions are entirely non-perturbative objects. But in the large transverse momentum 
region /Ca,fe± ^ ^qcd, "we can calculate the transverse momentum dependence, as shown in 
Sec. Ill, To compare to the results from the collinear factorization calculation, we let one 
of the transverse momenta ka± , kb± , be of order q±, and the others are much smaller. 
After integrating the delta function, we will obtain 

W'r" = ^ E F / TT {V.l-O^d - «) + - «)} . (81) 

W'aa - ^AA _ CisM V- f dx dz 



,2 



LT 



W^aa' + ^AA CisM ^ el f dx dz 



2 7r2Q3 

It is evident that the above results reproduce the differential cross sections we derived in 
the collinear factorization framework. 



VI. SUMMARY 

In this paper, we have calculated the naive-time-reversal-even but k±-odd TMD distribu- 
tions giT{x,k±), hiL{x,k±) at large transverse momentum, and they are related to a class 
of collinear twist-three matrix elements. We further studied the angular distribution of the 
lepton pair produced in the polarized Drell-Yan process for the single spin asymmetry Aut 
and double spin asymmetry Alt using the higher twist collinear approach. By comparing 
these results with those from the transverse momentum dependent approach, we found that 
they are consistent in the intermediated transverse momentum region. 

These calculations are not straightforward extensions of the previous calculations for the 
naive-time-reversal-odd TMD quark distributions [s?]. This is because, in the previous case, 
a pole contribution has to be taken in the final results, which will simplify the calculations. 
In this paper, we have to deal with more complicated kinemat ics, similar to next-to-leading 
order perturbative calculations for the gx structure function [42|. To carry out the calcu- 
lations, we have set up the twist expansion framework, and in particular, we have shown 
that the contributions from twist-three matrix elements will combine into the gauge invari- 
ant form. This shall encourage further developments in the higher-twist calculations. For 
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example, an extension to calculate the remaining TMD distributions hij,{x, k±) would be 
possible, though it will be more complicated because it is related the twist four collinear 
matrix element. The formalism we developed in this paper can also be extended to other 
semi-inclusive processes, such as the semi-inclusive deep inelastic scattering and back-to- 
back two hadron production in e'^e" annihilation processes. We will address these issues in 
future publications. 
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